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ABSTRACT 

The possibility that the strength of gravitational interactions might slowly increase with dis- 
tance, is explored by formulating a set of effective field equations, which incorporate the gravita- 
tional, vacuum-polarization induced, running of Newton's constant G. The resulting long distance 
(or large time) behaviour depends on only one adjustable parameter £, and the implications for the 
Robertson- Walker universe are calculated, predicting an accelerated power-law expansion at later 
times t ~ £ ~ l/H. 
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In the Standard Model of particle interactions, all gauge couplings are known to run with 
energy. Recent non-perturbative studies of quantum gravity have suggested that the gravitational 
coupling too may depend on a scale related to curvature, and therefore macroscopic in size. In this 
Letter, we investigate the effects of a running gravitational constant G at large distances. This scale 
dependence is assumed to be driven by gravitational vacuum polarization effects, which produce an 
anti-screening effect some distance away from the primary source, and therefore tend to increase 
the strength of the gravitational coupling. A power law running of G will be implemented via 
manifestly covariant nonlocal terms in the effective gravitational action and field equations. 

We start by assuming [1, 2] that for Newton's constant one has 

G{r) = G(0) [ 1 + c € (r/0 1/v + 0{(r/^) ] , (1) 

where the exponent v is generally related to the derivative of the beta function for pure gravity 
evaluated at the non-trivial ultraviolet fixed point. Recent studies have v~ l varying between 3.0 and 
1.7 [2, 3, 4, 5, 6]. These estimates rely on three different, and unrelated, nonperturbative approaches 
to quantum gravity, based on the lattice path integral formulation, the two plus epsilon expansion of 
continuum gravity, and a momentum slicing scheme combined with renormalization group methods 
in the vicinity of flat space, respectively. In all three approaches a non-vanishing, positive bare 
cosmological constant is required for the consistency of the renormalization group procedure. The 
mass scale m = £ _1 in Eq. (1) is supposed to determine the magnitude of quantum deviations from 
the classical theory. It seems natural to identify l/£ 2 with either some very large average spatial 
curvature scale, or perhaps more appropriately with the Hubble constant (as measured today) 
determining the macroscopic expansion rate of the universe, via the correspondence 

£ = 1/H , (2) 

in a system of units for which the speed of light equals one. A possible concrete scenario is one in 
which £ _1 = .ffoo = lim^oo H(t) = y/Tl\Ho with = ^, where A is the observed cosmological 
constant, and for which the horizon radius is H^-. 

As it stands, the formula for the running of G is coordinate dependent, and we therefore replace 
it with a manifestly covariant expression involving the covariant d'Alembertian operator 

□ = gTVuVv , (3) 



whose Green's function in d spatial dimensions is known to behave as 

5 \V > S(r-d(x,y\g)) ~ 



< x\^-\y > 5(r - d(x,y\g)) ^ > ( 4 ) 



where d is the minimum distance between points x and y in a background with metric g^v- We 
therefore write, in four dimensions, 



G -► G(D) = G(0) 



(5) 



One way of incorporating this is to replace the gravitational action 

1 



16vrG 

by 



dxy/gR (6) 



b / dx ^ ( x - CD fe) ^ + °^ u )~ 1,v ) ) fl 



The above prescription has in fact been used successfully to systematically incorporate the effects 
of radiative corrections in an effective action formalism [7, 8]. It should be noted that the coefficient 

in Eq. (1) is expected to be a calculable number of order one, not necessarily the same as the 
coefficient Co, as r and l/y/o are clearly rather different entities to begin with. One should recall 
here that in general the form of the covariant d'Alembertian operator □ depends on the specific 
tensor nature of the object it is acting on. 

The details of the incorporation of this modified G in the gravitational side of the Einstein 
equations are given elsewhere [10]. Here we shall describe instead its incorporation on the matter 
side of Einstein's equations, giving the effective field equations 

R»u -\g^R + Ag»» = 8nG (l + A(o)) , (8) 

where we have replaced G(r) by G(0)(1 + A(D)). These can be written in the form 

Rnu - \g^R + Ag^u = , (9) 

with T^ v = (1 + A(D)) defined as an effective, or gravitationally dressed, energy-momentum 
tensor. Just like the ordinary Einstein gravity case, in general T^ v might not be covariantly con- 
served a priori, V M T^ v / 0, but ultimately the consistency of the effective field equations demands 
that it be exactly conserved in consideration of the Bianchi identity satisfied by the Riemann tensor. 
The ensuing new covariant conservation law 

Vf^ = V [(l + A(U)) T^] = (10) 

can be then be viewed as a constraint on T^ v (or T^ v ) which, for example, in the specific case of 
a perfect fluid, will imply again a definite relationship between the density p(t), the pressure p(t) 
and the Robertson- Walker scale factor R(t), just as it does in the standard case. 



For simplicity we set the cosmological constant A to zero from now on and consider first the 
trace of the effective field equations 



R = 8vrG (1 + A(n)) T/ . (11) 

The advantage of this is that, initially, we need to consider the action of □ only on a scalar function, 
S(x) say, which is given by 

-^^Vs^S(z) . (12) 
V 9 

For the Robertson- Walker metric, 

ds 2 = _ dt 2 + R 2 {t) j_ilf_ +r 2 (^2 + sin 2^^| ? (13) 

□ acting on a scalar function of t only, is 



1 d 



R 3 (t) dt 

The energy-momentum tensor for a perfect fluid is given by 



(14) 



T^u = (p(t) + p{t)) Ulx u v + g^ P (t) . (15) 

We consider a pressureless fluid p(t) = and assume the density and scale factor are given by 
powers of t, as in the classical solution for the RW metric: p(t) = pot 13 , R(t) = Rot a . Then 

□«(W) = u n {-p{t)) 4"(-l) w+1 n {2 1 K R J J PotP- 2n . (16) 

v w; v y r(f + 1 - n)r(£±^±i - n) 

We may analytically continue the exponent to negative fractional n, [9] and obtain with n = 
— l/(2u), an expression for (1 + -A(D)) acting on the trace of T^, given by 

-(l + c< Q) V ^ po^ , (17) 

with 

c _ 4 -l/2i/ /_ 1 xl-l/2i/ ^2 + 1 Ml 2 J (18) 

Using the value of the scalar curvature for the Robertson- Walker metric in the k = case, 

R = 6(i? 2 (i) +R(t)R(tf) /R 2 {t) , (19) 

gives 

6a(2a-l) / fl , v 
^ " = " ( ! + c 5 (j) 1 Po* ■ (20) 
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For large t, when the correction term starts to take over, we see from the powers of t that 

(3 = -2-1/iz . (21) 

Next we will examine the full effective field equations (as opposed to just their trace part) of 
Eq. (8) with A = 0, 

Ry,v -\g^R = 8irG (1 + A(D)) . (22) 
Here the d'Alembertian operator 

□ = g^^pSv (23) 

acts on a second rank tensor, 

V u T a/ 3 = d v T a p - V\ v T\(i — Tp u T a \ = I va p 

V M (V„T a/3 ) = d^I ua p - - - T^I va \ , (24) 

and would thus seem to require the calculation of 1920 terms, of which fortunately many vanish by 
symmetry. Assuming that T^ u describes a pressureless perfect fluid, we obtain 



= r^( 2p(t)i ^ )2 ) 

= r 2 (l - kr 2 ) (DT^) rr 



(□T^) w = r 2 (l - kr 2 ) sin 2 6 (a T^) rr , (25) 

with the remaining components equal to zero. Note that a non-vanishing pressure contribution 
is generated in the effective field equations, even if one assumes initially a pressureless fluid. As 
before, repeated applications of the d'Alembertian □ to the above expressions leads to rapidly 
escalating complexity (for example, eighteen distinct terms are generated by D 2 for each of the above 
contributions), which can only be tamed by introducing some further simplifying assumptions. In 
the following we will therefore assume as before that k = 0, p(t) = p^t 13 ', and R(t) = R t a . We 
obtain 

(□ T^) tt = (6 a 2 - (3 2 - 3 a /? + /?) Po t^ 2 

(□T^) rr = 2R 2 t 2a a 2 Po t^ 2 , (26) 

which again shows that the tt and rr components get mixed by the action of the □ operator, and 
that a non-vanishing rr component gets generated, even though it was not originally present. 
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The geometric side of the gravitational field equations, the Einstein tensor, has the following 
components for the RW metric: 



G tt = 3R 2 (t)/R 2 (t) 

G rr = 1 ~l 2 (R 2 it) + 2 R(t) R(t 

G ee = r 2 (l-kr 2 )G rr 



G w = sin 2 6 Gee ■ (27) 

Then with k = and R(t) = Rot a , these will all behave like t~ 2 so in fact a solution can only be 
achieved at order provided the exponent (3 satisfies [3 = — 2 + 2n, or since n = —l/(2v), 

P = -2-1/v , (28) 

as was found previously from the trace equation, Eqs. (11) and (21). 

We must now determine a. By repeated application of □, for general n one can then write 

( an T IIU ) tt -► c tt (a,v) Po t~ 2 (29) 

and similarly for the rr component 

(O n T^) rr - c rr (a,v)R 2 t 2a Po t- 2 . (30) 

But remarkably (see also Eq. (25) ) one finds for the two coefficients the simple identity 

c rr (a,v) = \c tt (a,v) (31) 

as well as cee = r 2 c rr and c w = r 2 sin 9c rr . Then for large times, when the quantum correction 
starts to become important, the tt and rr field equations reduce to 

3a 2 r 2 = 8irGc tt (a,v) p t- 2 (32) 

and 

- a (3 a - 2)Rlt 2a ' 2 = 8nG c rr (a,v) R^t 2 * p t~ 2 (33) 

respectively. But the identity c rr = | c tt implies, simply from the consistency of the tt and rr 
effective field equations at large times, 

Crr(a,v) = i = 3a - 2 

c tt (a,u) 3 3a ' 1 ' 
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whose only possible solution finally gives the second sought-for result, namely 

a = \ . (35) 

We still need to check whether the above solution is consistent with covariant energy conserva- 
tion. With the assumed form for T^ v it is easy to check that energy conservation yields for the t 
component 

(V (□ T %„)) t - - ((3a + p + lMctt + Zacrr) potP+V"- 1 = (36) 

when evaluated for n = —l/2v, and zero for the remaining three spatial components. But from 
the solution for the matter density p(t) at large times one has f3 = — 2 — l/i/, so the above zero 
condition gives again c rr /c u = —(3a — 2)/3a, exactly the same relationship previously implied by 
the consistency of the tt and rr field equations. 

Let us emphasize that the values for a = 1/2 of Eq. (35) and (3 = — 2 — l/v of Eq. (28), 
determined from the effective field equations at large times, are found to be consistent with both 
the field equations and covariant energy conservation. More importantly, the above solution is also 
consistent with what was found previously by looking at the trace of the effective field equations, 
Eq. (11), which also implied the result (3 = — 2 — Eq. (21). 

The classical unmodified matter-dominated RW equations have solutions a = 2/3, (3 = —2, 
which mean that the scale factor behaves as 

R(t) ~ t a ~ t 2/3 (37) 

and the density as 

p(t) ~ t 13 ~ t~ 2 ~ {(R{t)y 3 . (38) 

This will also be the behaviour for our model at early times, but at later times, when the effect 
of the gravitational vacuum-polarization modification dominates, the behaviour is rather different: 
for the scale factor, we have 

R{t) ~ t a ~ i 1/2 (39) 

and for the density 

p{t) ~ IP ~ t" 2 " 1 ^ ~ (i?(t))- 2 ( 2+1 ^ . (40) 

Thus the density decreases significantly faster in time than the classical value (p(t) ~ t~ 2 ), a 
signature of an accelerating expansion at later times. 
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Within the Friedmann-Robertson- Walker (FRW) framework the gravitational vacuum polar- 
ization term behaves in many ways like a positive pressure term. The value a = 1/2 corresponds 
to uj = 1/3 in 

a = WT^) ' (41) 

(this follows from the consistency of the rr and tt equations in the general case) where we have taken 
the pressure and density to be related by p(t) = uj p(t), which is therefore characteristic of radiation. 
One can therefore visualize the gravitational vacuum polarization contribution as behaving like 
ordinary radiation, in the form of a dilute virtual graviton gas. It should be emphasized though 
that the relationship between density p(t) and scale factor R(t) is very different from the classical 
case. 
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